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The invariant mass of a jet is a benchmark variable describing the structure of jets at the LHC. 
We calculate the jet mass spectrum for Higgs plus one jet at the LHC at next-to-next-to-leading 
logarithmic (NNLL) order using a factorization formula. At this order, the cross section becomes 
sensitive to perturbation theory at the soft rn^ct/P^T^ scale. Our calculation is exclusive and uses the 
1-jettiness global event shape to implement a veto on additional jets. The dominant dependence on 
the jet veto is removed by normalizing the spectrum, leaving residual dependence from non-global 
logarithms depending on the ratio of the jet mass and jet veto variables. For our exclusive jet cross 
section these non-global logarithms are parametrically smaller than in the inclusive case, allowing 
us to obtain a complete NNLL result. Results for the dependence of the jet mass spectrum on the 
kinematics, jet algorithm, and jet size 7? are given. Using individual partonic channels we illustrate 
the difference between the jet mass spectra for quark and gluon jets. We also study the effect 
of hadronization and underlying event on the jet mass in Pyti-iia. To highlight the similarity of 
inclusive and exclusive jet mass spectra, a comparison to LHC data is presented. 
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I. INTRODUCTION 

There has been a rapidly expanding theoretical and 
experimental effort on techniques that exploit the sub- 
structure of jets (for a recent review see Ref. [l|). Jet 
substructure is of interest both for testing QCD and for 
identifying new physics. Much of the excitement in this 
field has been driven by the excellent performance of the 
ATLAS and CMS detectors, and the sophisticated jet 
measurements this has made possible at the LHC. Jet 
substructure measurements can for example be used to 
tag boosted heavy particles, whose decay products get 
coUimated into a fat jet, or to test and tune Monte Carlo 
programs. Most theoretical work has focused on design- 
ing these techniques and obscrvables with the help of 
Monte Carlo programs. At the same time, one would 
also like to know that these methods are under theoreti- 
cal control and build confidence that higher-order effects 
are not significant. (For some recent progress in this di- 
rection see e.g. Refs. 0-01 •) 

As our underlying hard process we consider pp ^ H +1 
jet with gluon fusion gg — > as the underlying Higgs 
production mechanism. This process is convenient as it 
provides a clean setup with a single quark or gluon jet 
in the final state via the three basic partonic channels 
gg — >■ Hg, gq — > Hq, and qq — >■ Hg. Of course, it is also 
important in its own right for Higgs measurements at the 
LHC, which rely on exclusive jet channels. 

Here we focus on one of the simplest jet substructures: 
the invariant mass of a jet. A successful calculation of 
this benchmark observable will instill confidence in our 
ability to carry out analogous calculations for other more 
complicated jet substructure observables. Such analyses 
require incorporating both a resummation of large loga- 
rithms al ln-'(mj/pj.^) where mj is the jet mass and pj. 
is the transverse momentum of the jet, as well as fixed- 
order perturbative corrections. This is made intricate by 



the dependence on multiple variables. There has been a 
lot of recent work on the calculation (resummation) of 
the jet invariant mass spectrum for jets with a realistic 
angular size 0, Isl-fisj which we will review in more detail 
below. Some of the key theoretical issues that must be 
addressed for the LHC case include: 

• Impact of summing large logarithms, \n{m^j/p'^^) 

• Soft radiation effects at the scale m?j/p^ 

• Impact of initial-state radiation 

• Color flow and hard process dependence 

• Dependence on kinematics including rapidity cuts 

• Jet algorithm and dependence on jet size R 

• Inclusive (> TV jets) vs. exclusive (= A'' jets) 

• Impact of non-global logarithms (NGLs) 

• Effect of hadronization on the spectrum 

• Effect of underlying event on the spectrum 

• Effect of pile-up on the spectrum 

• Utility of using groomed jets with trimming [l6j . 
filtering [l3|, or pruning [ll] 



We now elaborate on several of these items. For a jet 
with ^ 300 GeV, the jet mass peaks at mj ~ 50 GeV, 
leading to large logarithms of mj/pj.^ ^ 36. Therefore, 
a description of the peak region of the jet mass spectrum 
requires the all-order resummation of these logarithms. 
Soft radiation with momentum ~ m'^j/p^ is gener- 
ated by both initial and final-state particles and con- 
tributes at leading order in the power expansion to the jet 
mass. Since fixed-order corrections start to become rele- 
vant for resummation at next-to-next-to-leading logarith- 
mic (NNLL) order, a proper treatment of the soft scale 
~ m J /p'^ is crucial at this order [2, 19 . Numerically, 
the importance of these fixed-order soft corrections is also 
well known from recent work up to N'^LL for event shapes 
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in e+e" jets |15|. |22h24| . For processes with > 2 jets at 
hadron colliders there are multiple color structures, and 
the corresponding color flow must be taken into account 
starting at next-to- leading-logarithmic (NLL) order [2^ . 

The available freedom in defining a jet introduces a 
dependence of the jet mass spectrum on the choice of 
algorithm/clustering method and the jet size parameter 
R. There is also a choice of whether to use an inclusive 
or exclusive jet cross section, where the latter involves 
a veto on additional jets. The inclusive case has been 
studied at the LHC [26|, and inclusive calculations tend 
to focus on the anti-kT algorithm (Use of the anti- 
jet al gori thm avoids issues associated to clustering ef- 
fects |28l - l3ll |.) As wc will emphasize further below, a key 
difference between the inclusive and exclusive cases are 
the form of the non-global logarithms [H, [s^ that arise 
at 0{a1) beyond the Born cross section due to multiple 
restrictions on phase space. 

Let us summarize how the above issues have been stud- 
ied so far in the literature on jet mass calculations. The 
first calculations were carried out for event shapes in 
e^e^ — >■ jets using hemisphere jet masses. Here factoriza- 
tion theorems are well established and calculations exist 
up to N^LL [H [H, m In Refs. [1, Q a fac- 

torization formula for exclusive iV-jet cross sections at 
e"^e^ colliders was derived, where the angularity of a jet 
(which includes the jet mass as a special case) is mea- 
sured. This result only depends on the class of the jet 
algorithm (such as cone or kr-type), but suffers from 
non-global logarithms involving the jet veto and jet size 
R. The rcsummation of the jet mass in e+e~ — 2 jets 
with a jet veto was carried out at NLL in Rcf. , which 
includes a resummation of NGLs in the large- iVc approx- 
imation. This same process was considered in Ref. (Tl| . 
where the dominant R dependence of asymmetric thrust 
(which is related to jet mass) was obtained using a refac- 
torization of the soft function. In Ref. [l3 , this refactor- 
ization was verified at 0{a^) and the leading NGLs were 
obtained at this order. 

For jet mass calculations in pp collisions one considers 
jets with large transverse momentum, p;^, and with ra- 
pidities rj'^ away from the beam axis. Recently, several 
inclusive jet mass calculations have been carried out [isl - 
[T5| . In Ref. [ll], the jet mass was calculated using only a 
jet function. This ignores important contributions from 
wide-angle soft radiation, which couples together multi- 
ple hard partons, depends on the choice of jet algorithm, 
and contains NGLs. In Ref. the jet mass in — > 2 
jets and Z+1 jet were calculated at NLL, including a rc- 
summation of NGLs in the large- A'^c approximation. Al- 
though this is an inclusive calculation (no jet veto), one 
should also note that hard emissions giving rise to addi- 
tional jets are beyond the NLL order considered. In this 
case the dominant effect of the NGLs is on the peak of 
the jet mass distribution. Another inclusive calculation 
of the jet mass was carried out to obtain partial NNLL 
results in Ref. [l^, by expanding around the threshold 
limit. Here dynamical threshold enhancement (38l - l40| 



was used to argue that additional hard emissions are sup- 
pressed. Although NGLs were not resummed, their size 
was estimated, and found to mainly effect the peak region 
of the jet mass, as in Ref. [14 1. 

Our calculation at NNLL is for the exclusive jet mass 
spectrum, so it is useful to highlight differences with the 
inclusive case. At NLL, for a given partonic channel and 
fixed momenta of the hard partons, the two cases simply 
differ by a multiplicative factor, except for their respec- 
tive NGLs. In both cases the leading NGLs involve terms 
of the form 



-Pcut 



(1) 



for the cumulant jet mass spectrum integrated up to 
TO™*. For the inclusive jet mass spectrum, pcut is a 
hard scale ~ p^ and the NGLs are therefore large loga- 
rithms that are parametrically of the same size as other 
a* In-* (toj/pj.'^) terms, and are thus part of the NLL re- 
sult. Hence, in this case a complete resummation at 
NLL (or beyond) requires the NGLs to be resummed 
to all orders, which practically is currently only possi- 
ble in the large- A^c approximation. In contrast, in the 
exclusive case pcut is an adjustable parameter and is re- 
lated to the jet veto. Here we have both TOj ^ p;^^ and 
Pcut ^ P't^i so the logarithms in Eq. ^ are smaller than 
in the inclusive case. In particular, for fixed Pcut there is 
a point in the uij spectrum where the NGLs vanish, and 
there is a region about this point where the NGLs are 
not large logarithms. An estimate for the size of this re- 
gion can be obtained from the series of three NGL terms 
(log-squared, log, and non-log) that are known for the 
hemisphere jet masses [4l|, . When all the terms in 
this series are of similar magnitude the logarithmic en- 
hancement is not dominant, and the NGLs do not need 
to be resummed. This occurs for 1/8 < TOj"*^/p^^^ < 8. 
We will numerically explore the size of this region in our 
exclusive jet mass calculation, and demonstrate that the 
region is large enough that we may consider the non- 
global logarithms to not be large. This can be contrasted 
with Fig. 3 of Ref. [l3|, which shows that the presence of 
an unmeasured region of phase space makes large NGLs 
unavoidable in the inclusive case [13, [Hi . 

It should also be noted that although exclusive jet cross 
sections are not necessary for jet mass spectra, they are 
important in their own right because many Higgs and 
new physics searches categorize the data by the number 
of jets to improve their sensitivity. For example, the im- 
portance of the Higgs + 1 jet channel in if —)■ tt and 
H WW* was pointed out in Refs. [4l,[43|. Recently a 
NLL rcsummation of jet veto logarithms was carried out 
in the context of Higgs plus jets in Ref. [i^ . 

Our calculation of the jet mass is centered on using 
the A^-jettiness global event shape [1^ to define jets, in- 
stead of a more traditional jet algorithm. For an event 
with A^ jets, A^-jettiness assigns all particles to A^ -I- 2 
regions, corresponding to the A^ jets and two beams. We 
calculate the cross section for pp —)■ -I- 1 jet at NNLL, 
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fully differential in the contributions of each region to 1- 
jettiness. For the jet region, this contribution yields the 
jet invariant mass. The contribution from the remaining 
two beam regions are used to implement the jet veto. In 
each of these variables there is a series of large double 
logarithms that must be summed. 

An advantage of using iV-jettincss is that the jet veto 
is made through a jet mass-type variable, rather than a 
Pt variable. Therefore, the structure of the perturbation 
theory, which is simultaneously differential in these two 
kinematic variables, is simpler. In particular, there is 
a QCD factorization formula for this cross section [i^, 
|47| . obtained by making use of Soft-Collinear Effective 
Theory (SCET) [l^Iii-^. For the experimentally more 
realistic case of measuring mj with a pt veto variable one 
must simultaneously deal with a thrust-like resummation 
and a pr-type resummation. 

Returning to our list of theoretical issues from the be- 
ginning, the use of iV-jettiness allows us to carry out 
the summation of large logarithms at NNLL while prop- 
erly accounting for soft radiation effects and initial-state 
radiation. We also use it to calculate the dependence 
of the jet mass spectrum on the jet kinematics, the jet 
size, and the definition of the jet region. Results are 
shown for individual partonic channels, gg — >■ Hg and 
gq —J- Hq, illustrating the differences between quark and 
gluon jets, as well as the full pp H + 1 jet process 
from the Higgs coupling through a top quark loop. To 
investigate the differences between exclusive and inclu- 
sive jet mass measurements we compare our results with 
Pythia and also to ATLAS jet mass data [l^ . We also 
analytically explore the effect of NGLs on the jet mass 
spectrum, and the effect of hadronization and underlying 
event with Pythia [51|, [12] . 

Thus, we address all items in the list of issues except 
for the last two, for which some brief comments are in 
order. Methods for removing pile-up contributions to 
jet observables have been discussed in e.g. Refs. [ssl. [s^. 
and direct pile-up calculations are beyond the scope of 
our work. Finally, it is known that grooming jets has a 
large impact on their soft radiation and causes significant 
changes to the jet mass spectrum. We do not attempt to 
analytically control the effects of jet-grooming methods 
here. 

In calculating the jet mass we consider both absolute 
and normalized spectra. Normalizing the jet mass spec- 
trum reduces the perturbative uncertainty, and turns out 
to remove the dominant dependence on the jet veto vari- 
able. In particular, the jet veto dependence cancels up 
to NLL if we consider a particular partonic channel and 
fixed jet kinematics. We will show that this cancellation 
remains effective when summing over partonic channels 
and integrating over a range of jet momenta. 

In Sec. ini we discuss the kinematics and several jet 
definitions based on iV-jettiness, exploring their features. 
The technical details of our calculation are presented in 
Sec. mil Here we discuss the factorization formula for 
the cross section, the rcfactorization of the soft func- 



tion, non-global logarithms, and the choice of running 
scales. Sec. IIVI and Sec. |V] contain our numerical results 
for the individual partonic channels and ioi pp ^ H + 1 
jet, showing the dependence of the jet mass spectrum 
on the jet veto cut, the order in perturbation theory, 
the jet kinematics, the jet definition, the jet area, on 
gluon versus quark jets, and on NGLs. Using PythiaS, 
in Sec. IVII wc analyze the hard process dependence for 
gluon jets, compare inclusive versus exclusive jet mass 
spectra, study the dependence on classic jet algorithms, 
and look at the impact of hadronization and underlying 
event. We also compare our NNLL exclusive jet results 
with Pythia for the same jet definition and kinemat- 
ics, and compare them with inclusive jets from the LHC 
data. We conclude in Sec. IVIII Detailed ingredients for 
the NNLL cross section are summarized in appendices. 

II. KINEMATICS AND JET DEFINITIONS 

We describe the process pp ^ H + 1 jet using the 
transverse momentum p'^ of the jet, the pseudorapidity 
rjj of the jet, and the rapidity Y of the hard collision 
relative to the CM frame of the colliding protons. The 
1-jettiness event shape is defined as [i^ 

^ . r 2(?,7 • pk 2qa ■ Pk 2qb ■ Pk \ 

k 

where a, b denote the two beams and J the jet, the qi are 
massless reference momenta and the Qi are normalization 
factors. For the reference momenta we take 

qj = Ejil,rij), g^,, = a;a,6^(l,±z) . (3) 

The jet energy Ej and jet direction nj can be predeter- 
mined with a suitable jet algorithm. The jet algorithm 
dependence this induces on 7i is power suppressed [i^, 
and we will use anti-kT.-'^ The unit vector z points along 
the beam axis, and the momentum fractions Xa and Xb 
arc fixed in terms of the total invariant mass Q and ra- 
pidity Y , 

XaXbEl^ = = {q.j + qnY , 

ln^=2y.ln%-^f/^^, (4) 
Xb [l,z) ■ (q,j + qH) 

where denotes the momentum of the Higgs. For later 
convenience we also introduce the notation 

Sij 2gfj • qj . (5) 



If Qj = 2-Ej then an equally good choice would be to minimize 
Tl with respect to the axis nj. A fast algorithm to carry out 
this minimization has been devised in Ref. [55j . using a slightly 
different A'^-jettiness measure than the ones we use here. 
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(a) Geometric E and pT measures with p = 0.5. 
FIG. 1: Comparison of the jet regions for different jet meat 

The minimum in Eq. ^ divides the total phase space 
into 3 regions, one for each beam and one for the jet. We 
denote their contributions to 71 as 7^ and 71 for the two 
beam regions, and 7j for the jet region, so 

Ti=Tj + Ta+Tb. (6) 

The contribution of the jet, 7j, is directly related to the 
jet's invariant mass mj 

= -Pj = ("./ ■ Pj)inj ■ pj) - p]±_ 

= 2qj-pj [1 + 0(A2)] 

= QjTj[l + OiX')], (7) 

where pj is the full jet momentum defined by summing 
all particles in the Tj-region, n^j = (l,nj) and fij = 
(l,—nj) are defined by the predetermined jet direction 
nj, and the power counting parameter A scales as ^ 
Tj/Ej ~ rrij/Ej. In the second line of Eq. ^ wc used 
the fact that nj and the exact direction of the iV-jettiness 
jet, Pj, differ by very little, such that pj±/{nj -pj) ^ A^. 
The difference between these two jet directions affects 
the jet boundary, which changes the contribution of soft 
radiation to the jet pT, but only by a small amount ^ A". 
We also used that the large jet momentum fij -pj ~ nj ■ 
qj[l + C(A2)]. For a jet with - 300 GeV these 0{X^) 
power corrections are 1/36 ~ 3% in the peak region, and 
hence negligible relative to the perturbative uncertainties 
at NNLL. Investigating the jet mass spectra for the exact 
TO J — -P/ using TOj — QjTj in Pythia, wc also find 
that they arc indistinguishable. 

The details of the beam and jet regions selected by the 
minimum condition in Eq. ([2|) depend on the normaliza- 
tion factors Qi. Since their values affect which particles 
are grouped into the beam and jet regions, they con- 
stitute a jet measure. They also impact the geometric 
shape of the jet area. Differences between measures are 
therefore similar to the different choices for jet-algorithms 
(anti-kr, Cambridge-Aachen, cone, etc.). We will con- 
sider a variety of choices: 




(b) Anti-ky and geometric R for R = 1. 
3 at different rj and (p. The marks the jet direction nj. 

• invariant-mass measure: 

Qj = Qa = Qb = Q (8) 

• geometric pT measure: 

Qj = 2p I^tI = 2pEj/ coshryj (9) 

Qa^b — ^a,6-£^cm — Q 

• geometric E measure: 

Qj^2pEj (10) 

Qa,b ' — ^a,bEcm — ^ Q 

• geometric R measm^e: 

Qj^2p{R,7]j)Ej (11) 

Qa^b '^a^bEcin C Q 

where p{R,rij) fixes the area of the jet in {ri,(j))- 
spacc to be ttR^. 

In all cases p is a dimcnsionless parameter that allows 
one to change the size of the jet region. In the geometric 
R case p is fixed in terms of the jet radius parameter R.^ 
The choice of Qa.b hi the geometric measures removes the 
dependence in q^/Qa and q'^ /Qb on the total rapidity Y . 
This is useful in the presence of missing energy, which 
prohibits the measurement of the boost Y of the partonic 
center-of-mass frame. Since for the geometric measures 
Qj ^ Ej, they are all insensitive to the total jet energy. 
For the geometric px case wc have explicitly 

=PkT (2^^^^ coshAyifc-2cos A^ifc) (12) 

QiT ^ PkT ' 



^ For the multijet case we would use the same p{R,r]j) for each 
jet that is determined when they do not overlap. 
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FIG. 2: Numerical results for p{R,rij) in the geometric 7? measure. Left: Dependence on R for rij = 0, which is ~ B? as 
expected. Right: Dependence on r]j for R — 1. To solve for p we use a fit (solid line) to the true rij dependence (dots). 



where Ay^fc =yi-yi,, A(j).ik = (t)i-(pk are the differences 
in rapidity and azimuthal angle between the direction of 



jet i and particle k, and Tnj.rp 



pIt 



for a particle 



of mass ?n. For masslcss particles wc thus get 
2qi ■ Pk 



QiT 



PkT (2 cosh Ayik - 2 cos A^^^) 
VUT [(Ay)2 + (A0,,)2] . 



(13) 



The jet regions for geometric pT and E are roughly circu- 
lar, as shown in Fig. 1(a) They become smaller at large 



rapidities for geometric pt, while they stay of comparable 
size for geometric E. 

For geometric i?, numerical results for the parameter 
p(i?, 77,/) as function of R and 77,/ are shown in Fig.[2l The 
left panel shows that the dependence on the jet radius R 
is approximately p ^ R? , as expected. The right panel 
illustrates the dependence on r]j for fixed i? = 1 , showing 
that p approaches a constant for large 777, i.e. when the 
jet becomes close to the beam. When using geometric 
R in our results below, we use for convenience a fit of 
the 77j dependence for fixed value of R. For example, for 
R = 0.5, 1, 1.2 we have for |77,7| < 2 



Note that for R = 0.5 the parameter p increases rather 
than decreases with r/j. A comparison of the jet regions 
for geometric R with anti-kr jets is shown in Fig. 1(b) 



Although their areas are chosen to be the same, the geo- 
metric R jets are not perfectly circular and have an "off- 
set" between the jet direction and the center of the jet 



region. The former (latter) effect decreases (increases) 
with |77j|. For a smaller jet radius of i? = 0.5 the geomet- 
ric R jets become more circular also at central rapidities 
and very close to anti-ky jets. In Ref. [5^ a modifica- 
tion of iV-jettiness was introduced that matches anti-ky 
closely for any R. However, this definition reintroduces a 
region of phase space that belongs neither to the jet nor 
the beams, making it more complicated for calculations. 

III. CALCULATION 

A. Factorization Formula 

We start by rewriting the phase space integrals for the 
hard kinematics in terms of the rapidity 777 and trans- 
verse momentum of the jet and the total rapidity Y , 



dxb f d^QH f d^ 



J Eh J Ej 
= /dr/jdp^dr ' 



i27r)^S\qa + qb 



qj 



1 p-^ 



2tt Q2 



1) 



(15) 



The variables were defined in Sec. |TI1 and we used az- 
imuthal symmetry and the relations 



piR = 


0.5, 77j) 


= 0.164 + 


0.03777^^ - 


- 0.00977) 


+ 0.000877^ , 




p{R 




= 0.834- 


0.23377} - 


h 0.07777;^ 


- O.OO8776 , 


Pt 


p{R = 


1.2,77,7) 


= 1.272- 


0.37677} - 


h 0.10177;^ 


- O.OIO77*}. 

(14) 


Q 



H 



2Q cosh(7/j - Y) ' 

h(77,/ -Y) + ^p-p cosh^(77j -Y) + mj^ 



(16) 



p'^ cos 



Most of our plots will be normalized and for fixed values 
of 77,/, p'^, and y, in which case the phase space factor in 
Eq. pS]) drops out. 

Our calculation relies on the iV-jettiness factorization 
formula in Ref. 47| , which we here specialize to the case 
of 1-jcttiness: 



d^a^+ij' _ pJ, 



-YT'^Hf^{{qf},p) dtaBf^^{ta,Xa,n) dtbB^^{tb,Xb,p) 



d77j dp^ dY dTa dTb dTj 47r£;2„^ (Q2 + ^2 

x/d.,j«,,(.„M)^.(r.-^,r.-^,r.-g,{|},^ 



p . (17) 
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channel 




Kb 




gg Hg 


9 


9 


g 


gq Hq 


9 


q 


q 


qg Hq 


q 


9 


q 


gq -i- Hq 


9 


q 


q 


qg Hq 


q 


9 


q 


qq->- Hg 


q 


q 


g 


qq Hg 


q 


q 


g 



TABLE I: Values of n for the different partonic channels. 



The A^-jcttincss variables Ta, Tb, and Tj were defined in 
Sec-ini The hard function 77^ contains the short-distance 
matrix element for producing a Higgs plus a jet, the beam 
functions and B^,^ describe the coUinear initial-state 
radiation and contain the PDFs, the jet function J^j 
characterizes the collinear final-state radiation, and the 
soft function 5^ describes soft radiation effects.^ The 
sum over k = {Ka, Kb, kj} runs over the possible flavors 
Ki € {g, It, M, d, . . . } of the two incoming and one outgoing 
parton. The possible combinations, corresponding to the 
various partonic channels, are listed in Table HI 

The power of factorization is that it allows one to eval- 
uate the various fixed-order pieces at their natural scales, 
where they contain no large logarithms. We then use the 
RG evolution of each of these functions to evolve them to 
a common scale fj., resumming the logarithms of m^j/pjP' 
and QiTi/px^. This evolution is implicit in Eq. (fTT)) . by 
writing all functions as evaluated at the common scale 
fi. The factorization formula with all evolution factors 
written out explicitly is given in Eq. (|28p below. Our 
choice of scales is discussed in Sec. IIII CI Power correc- 
tions to Eq. (ini) arise from so-called nonsingular correc- 
tions, which are suppressed by a relative 0(rn?j / Q^) in 
this differential cross, and are not considered here. 

The cross section in Eq. p7|) is differential in the 1- 
jettiness contributions from the jet and the beams 7j, 
Ta, and Tb- As we will see, the shape of the jet mass 
spectrum is independent of the jet veto for a reasonable 
range ofTa.b values. For simplicity we impose a common 
cut Ta.b < T^"'. We also convert Tj to the jet mass mj 
using Eq. Q. and so consider 



T" 



dTa [ dTb 
Jo 



dTj 



dV 



dTa dTb dTj 



(18) 



The differential jet mass cross section, da/dm j, is ob- 
tained by taking the numerical derivative of this cumu- 
lant cross section. We define the normalized jet mass 





matching 


7x Tcusp P PDF 


LL 


0-loop 


1-loop 1-loop NLO 


NLL 


0-loop 


1-loop 2-loop 2-loop NLO 


NNLL 


1-loop 


2-loop 3-loop 3-loop NLO 



TABLE II: Perturbative ingredients at different orders in re- 
summed perturbation theory. 



spectrum over the range [0,™™*] as da/dmj, so 



da 
dm J 



1 



da(r="*) 



r™*) dm J 



(19) 



The ingredients in the resummed cross section are 
needed at different orders in perturbation theory, as sum- 
marized in Table im where the columns correspond to the 
fixed-order matching, non-cusp anomalous dimension , 
cusp anomalous dimension Fcuspi the /3 function, and the 
PDFs. All ingredients necessary for a NNLL resumma- 
tion of the global logarithms arc known and are collected 
in App.El The one-loop hard function for the three basic 
processes gg — Hg, gq — >■ Hq, and qq — >■ Hg via gluon 
fusion (in the large mj limit) are obtained from the one- 
loop helicity amplitudes calculated in Ref. [11] following 
the procedure in Ref. [13 . The one- loop quark and gluon 
jet function were calculated in Refs. |58l- [6(][ . the one-loop 
quark and gluon beam functions in Refs. |6ll - [63 |. and the 
one- loop soft function in Ref. [13 . We also require the 
cusp anomalous dimension to three loops (65l . l66j . and 
the non-cusp anomalous dimensions to two loops, which 
are known from Refs. [fiS [sl [gMII- 

There is some freedom in how to treat products of the 
fixed-order corrections in Eq. (|17p . specifically the higher- 
order cross terms that are generated, such as the one-loop 
correction to H times the one-loop correction to J, which 
we denote H'^^W-^\ The scries for the individual objects 
are fairly convergent, except for the hard function whose 
one-loop correction is known to be rather large. ^ We 
therefore expand the convolutions of the fixed-order B, 
J, and S functions order by order in as to the order 
needed, but keep the hard function i?^"^ + H^^^^ as an 
overall multiplicative factor. We anticipate that the cross 
terms involving if^^-* will account for a large portion of 
the NNLO cross section. Schematically, this means that 
the fixed-order components of our cross section take the 
form 

(i7(o)+ ^(D) [(5(0)^(0) j(o)) ^ ^(0) (20) 

+ (S(i)s(o) j(o)) ® ^(0) + (sWfid) j(o)) ® ^(0) 
+ (B^SW J(i)) ® + (S(°)S(°) J(0)) ® 5(1) 



3 Note that we do not call Eq. JTtJ a factorization theorem since 
the decoupling of Glauber gluons for hadron collider processes 
with a specific number of jets has not been proven. 



* We find that the improvement of using a complex scale fij{ for 
the hard function is only marginal, because some logarithms are 
minimized by a negative fi^ and other by a positive 
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The difference between expanding the hard function or 
treating it as multiphcative is within our pcrturbativc 
uncertainty, being a ^ 20% effect for the unnormahzcd 
mj spectrum, and only a ^ 2% effect for the norniahzed 
mj spectrum. 



B. Refactorization of the Soft Function 

For a process with one or more jets there are multiple 
directions for collincar radiation and various kinematic 
variables so a few additional hierarchies become possi- 
ble. The factorization formula assumes that there are no 
additional strong hierarchies beyond the collinearity of 
the jet TOj ^ Pt^, and the absence of additional central 
jets away from the beam directions, namely QaTa ^ P^^ 
and QbTb p^^- Physically, this corresponds to the 
following four assumptions 

1) QiTi ^ QjTj commensurate mj and jet veto 

2) ^ ji'' ^ ^ ji'' well separated jet and beams 
EiEj EiEk 

3) Ei ^ Ej jet and beam-jets of similar energy 



4) ^Q, 



jet and beam regions of similar size 



Assumption 1) ensures that we are in the region where 
NGLs arc not large logarithms. Assumption 2) implies 
that the jet is not too close to the beam direction, and 
avoids having large angular logarithms, which would re- 
quire an additional "ninja summation" Q. 

Three combinations of these four assumptions are nec- 
essary to avoid introducing additional large logarithms 
that are not summed by the renormalization group evo- 
lution of terms in the factorization formula, namely 



Sik 



Qi. 



(21) 



The first implies that the logarithms in the hard function 
can be minimized with a common scale /i, and all three 
combine to imply that a common scale also minimizes 
all logarithms in the soft function. One combination of 
assumptions, Ei/Qi ^ Ej/Qj, does not appear explicitly 
in arguments of functions in the factorized cross section, 
and hence does not show up in logarithms for the lead- 
ing power result. However, it is in general necessary as 
part of the derivation of Eq. (fT7|) to ensure that certain 
neglected terms are power suppressed. 

An important consideration in carrying out the sum- 
mation of large logarithms is the order in Us and loga- 
rithms at which violations of Eq. (|2ip first become ap- 
parent. For the soft function the first terms that appear 
for the various logarithms arc 



a, ln2(- 



a. In 



at W 



where we integrate the soft function over % up to the 
cumulant variable T^. The first of these is part of the LL 
series. The second is an angular logarithm. It is part of 
the NLL series if it counts as a large logarithm. Other- 
wise, it is part of the ^ as fixed-order terms that start to 
contribute at NNLL. The third is a NGL. It is part of the 
NLL series if it is a large logarithm. Otherwise it is part 
of the ~ a1 fixed-order terms that start to contribute 
at N"^LL. Therefore, there is a nontrivial constraint on 
the choice of scales ^ in the soft function. The scales 
must be chosen to minimize the first type of logarithm in 
Eq. (22) without inducing terms of the form of the second 
and third types already at LL order. In particular, this 
implies that a poor scale choice could introduce unphys- 
ical angular logarithms or NGLs into the LL series. For 
our choice of kinematics and Qi the second type of angu- 
lar logarithm in Eq. (|22p is never large. However, since 
we are exploring a spectrum in m?j = QjTj the third 
term in Eq. (|22p will grow as the parameters are varied. 
To surmount this problem requires a refactorization of 
the soft function which we will consider below. 

For the hard function the series of leading double log- 
arithms involves terms of the form 



a. In 



V s, o J 



an 



Sik 



(23) 



(22) 



For the choice of jet kinematics explored in this paper we 
will always satisfy the assumption s^j ^ s^fc, so there is 
no additional constraint on the scale associated with the 
hard function. 

The hierarchy between 7j and 7"'^"* leads to unphys- 
ical large logarithms if a single scale /xg is used for the 
initial conditions for the soft function evolution. Here we 
address how these can be removed by a refactorization of 
the soft function, with corrections from the true hig her 
order non-global logarithms (see Refs. [1, [H, [3 l37l | for 
earlier refactorization discussions). 

In general, the all-order soft function has the form 

^({fcJ,{C}>M) 

= n S^[h,M'l}. M) + 5NGL({fc,;}, t^) , (24) 

i 

where = q'^/Qi- Here S'ngl contains all non-global 
terms, and hence has an intrinsic dependence on the ra- 
tios ki/kj. At NLO there is only one soft gluon emitted, 
which can contribute to only one of the 7i at a time. Thus 
the NLO soft function factorizes, and S'ngl ^ C(«s)- 
Truncating to 0{as) there is still some freedom in the 
definition of the Si. Whereas the terms with explicit 
ki dependence in S{{ki},fi) clearly belong to Si{ki,ii), 
the pure delta function terms 6{kj)S{ka)S{kh) can in 
principle be split in multiple ways between the various 
Si{ki,iJ,). We choose to split these terms evenly, as de- 
tailed in App. IA41 and we introduce an additional pa- 
rameter r in the scale variation to estimate uncertainty 
from this freedom as discussed further below and in detail 
in Sec. ImCl 
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Due to the consistency of the factorization formula, 
the evolution of the soft function factorizes exactly to all 
orders in perturbation theory, 



(25) 



(Here we used the fact that the beam and jet functions 
have the same evolution [g^.) This involves the factor- 
ization of the evolution of the hard function H = CC^, 
which follows from the form of the anomalous dimension 
for C fzllTi, 



7c (m) = rcusp[as(/i)] 



V MO 



^T,.T,ln(-f^-iO 



+ 7c[«.(a^)]. (26) 



The sum on i and j runs over the colored partons partic- 
ipating in the short-distance interaction and denotes 
the corresponding color charge matrix. (For pp ^ H + I 
jet the color space is still trivial, so these color matrices 
are just numbers.) To associate the In fj. terms to individ- 
ual partons we introduced a dummy variable /io and used 
color conservation. It is not a priori clear how to associate 
the remaining terms within the X]i<j to each Us^, and 
we choose to split each term evenly between i and j. The 
explicit expression for the factorized hard function evolu- 
tion that we employ is given in App. lA 51 Other potential 



choices of splitting up these terms are again probed by 
the scale parameter r, which is discussed in more detail 
around Eq. (|5T|) . and the corresponding uncertainty is 
found to be small. The two-loop non-cusp anomalous di- 
mension has the structure 7c(Q!s) = nqjq + rig^g, where 
rig and Uq arc the number of gluon and (anti)quark legs, 
so it naturally factors. 

The factorization of the evolution and fixed-order soft 
function in Eqs. ([24|) and ((25|) suggests that we can eval- 
uate the piece of the soft function corresponding to % at 
a scale , 

S{{ki},^Ji) = n / '^KUsAh - kl^^,^,s.)S,{{k'^},fls^) . 

i 

(27) 

This factorization does not hold for all the terms at or- 
der a^, since there are diagrams that contribute to mul- 
tiple 7i, leading to non-global logarithms of the form 
a^ln^(fcf/fc|) appearing in S'ngl in Eq. p4|) . We dis- 
cuss in Sec. IIII Dl how we estimate the size of these NGL 
contributions in the jet mass spectrum. 

In our implementation we find it simplest to run the 
hard, jet, and beam functions, rather than the soft func- 
tion, as summarized in Fig. [31 though the final results 
are completely independent of this choice. Since the cut 
on both beams is the same, they have a common fiB and 
(J-Sb ■ We summarize the work in this section by pre- 
senting the factorization formula valid at NNLL which 
includes the evolution factors and refactorization of 5*, 



dr]j dp^ dV dmj At:E^^{Q'^ + m\ ) 

X / dta At'^ Uj^^ (ta - t'^, HSb > tJ-B)B^^ (t^, Zq, ^b) / dtfe At[ Uj (tb - 4, ^1SB : IJ'B)Bk, {t'b, Xb.Hs) 



All necessary perturbative results for H^, J^,, Si, and 
the Ui are collected in App. [A] 



C. Choice of Running Scales 

The factorization formula in Eq. sums the large 
logarithms of QiTi/pj^^ from the cuts on the beams and 
of QjTj Ipx^ = rrij/p'r^'^ from the jet mass measurement 



by carrying out perturbation theory for the hard, beam, 
jet, and soft functions at their natural scales 

Mff-PT7 fJ-t-V QiTi , ns.-QiTi/pT, (29) 

and then running them to an arbitrary common scale. 
Here i = a,b, J, since the situation for the beams and the 
jet are fully analogous. As discussed above we must use 
the combination QiH to ensure we have the correct lead- 
ing logarithms. The remaining dependence of the cross 
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FIG. 3: Illustration of the different fixed-order scales appear- 
ing in the factorized cross section and our evolution strategy. 



section on Qi occurs due to its impact on the boundaries 
between the jet and beam regions. This is encoded in the 
fixed-order terms in the soft function that do not involve 
large logarithms. 

If any 7i becomes very small, ^ Aqcd, the nonper- 
turbative corrections to the soft function become impor- 
tant and the scales are i^h — P'tj l^i ^ ^/Aqcd P'tj ^^id 



f^Si ^ Aqcd. On the other hand, for large % ^ the 
resummation of Qili/pj'^ is not important and is turned 
off by having the scales merge, = fJ-i = fJ-Si — Pt- 

To connect the various regions where the resummation 
must be handled differently, we use 7i-dependent scales, 
which arc known as profile functions (23l . iTSj . A transi- 
tion between these three regions is given by the following 
running scales, adopted from Ref. [64| . 



fJ-H = fJ-: 



(30) 



/Xi(T) = 1 + e,; e'(T3 - r) ( 1 ) \/ H Idrunir, /i) , 



T3 

MS, (t) = 1 + es. ^'(^3 Mrun(T, fl) , 

where r = QiTi/Q^- The expression for /irun can be 
found in App. [Bl along with the choice of parameters for 
the central values, and details on the variations of /i, e^, 
and used to estimate the perturbative uncertainties 
in our predictions. 

To estimate the additional perturbative uncertainty as- 
sociated with the refactorization of the soft function in 
Sec. IIIIB| we reintroduce correlations between the soft 
scales using a parameter r satisfying < r < 1, 



('•) 

ln/2s^^ 



In fJ-SB 



(r) 



Tj In i-J-Sj 



T2 



(31) 



Here = Cp for i = q and i = q, and = Ca for 
i = g. For r we have the original uncorrelated 
soft scales. By increasing r the scales move towards the 
"color average" value fls- At r = 1 they are all equal 
to this average soft scale, so the refactorization is turned 
off (which as explained earlier causes unphysical NGLs in 
the LL scries). To estimate the size of the freedom in the 



refactorization we take r = 0.2 as our default choice and 
include r = and r = 0.4 as separate scale variations in 
our uncertainty estimate. 



D. Non-Global Logarithms 

If the NGLs are not large logarithms then they enter 
beyond NNLL order, and should be of comparable size to 
other higher-order perturbative terms. This is obviously 
only possible for some range of mj/(pyT'^"'), which de- 
termines where our result is valid at NNLL order. To 
determine this range we include the leading O(a^) NGL 
into our resummed calculation and compare the results 
with and without this term for various parameter choices. 
In the factorized cross section all NGLs enter through the 
soft function S. For simplicity we restrict this analysis of 
the size of non-global logarithms to the gg — > Hg chan- 
nel, as the results for other channels arc similar. 

The leading NGL in the cumulant soft function is 



^NGL({fca,Ms)=n(/''' 



dk,)S^GL{{k.},fis) (32) 



(2^)2 



i<j J 



Here dj is a geometric factor that depends on the bound- 
aries of the jet and beam regions. Note the absence of 
explicit /i-dependence in the NGLs. These expressions 
follow from the known result for e~^e~ — > 2 




^,|4^, by replacing the color factor CfCa — >■ 
C^. Unlike the global logarithms this contribution does 
not factor, so we assign it a common soft scale which we 
take to be fls given in Eq. (|3T|) . 

For the purpose of our numerical analysis we take 
Gij = TT^/S, which is the result for a hemisphere. This 
may be thought of as reasonable estimate and in reality 
the values may differ by about 15% to 30% [l3]. Con- 
verting the cumulant space result in Eq. ([32|) into a full 
distribution yields 



«|(M5)^^ 

(2^)2 3 



(33) 



where the £„ denote standard plus distributions as de- 
fined in Eq. (jASp . Note that the /i' dependence cancels 
out explicitly between the terms, so the choice of this 
scale is arbitrary and irrelevant. It is introduced for cod- 
ing purposes, since it is convenient to have the same type 
of Cn distributions as in the non-NGL part of the soft 
function. When the NGLs are included in this manner, 
via the soft function in the factorization, one automati- 
cally resums an infinite series of global logarithms that 
multiply the NGL. In particular, this includes terms that 
arc schematically [a^ ^'0-^][J2ki'^s In^)'^] where the first In^ 
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is non-global and the second In is a large global loga- 
rithm. The all-order structure of this series of terms is 
correctly predicted by the factorization. 

For our analysis we will mostly be interested in the 
normalized spectrum in Eq. p9[) . Here in the numerator 
the two jet veto variables are in cumulant space and mj 
is in distribution space, while in the denominator all the 
variables are in cumulant space. This result has two types 
of NGLs 



i) 
ii) 



a^(/is)ln 



'J 



V 



(34) 



For the denominator the relevant form of the NGL log- 
arithms is as in Eq. (f5^ . yielding the terms i). For the 
numerator the form of the NGL is as in ii). The presence 
of two types of NGLs in the normalized spectrum implies 
a somewhat different dependence than for the unnormal- 
izcd cross section. The effect of NGLs in these two cases 
arc analyzed in detail in Sec. IIV CI There we will show 
that there is indeed a fairly large range of to™' values 
where the NGL terms in the exclusive jet cross section 
are not large logarithms. 



IV. RESULTS FOR GLUON AND QUARK JETS 

In this section we focus on the individual quark and 
gluon channels, leaving results for pp ^ H + I jet to be 
discussed in Sec. |V] below. We first study the theoret- 
ical predictions for the toj spectrum with and without 
normalization, and show that normalizing substantially 
reduces the perturbative uncertainty. We also study the 
order-by-order convergence of this differential cross sec- 
tion, and the size of various contributions to the pertur- 
bative uncertainty bands. Next, the dependence on the 
jet veto T"'^"' is studied. Finally, we investigate the size 
of non-global logarithms as a function of toj and T'^^*. 



A. Default Parameter Choices 

Unless indicated otherwise we use the following default 
parameter choices for all plots in Secs. lIVilyl and I VII For 
the Higgs mass we take mn = 125 GeV |74l.[75l|. and for 

= 7TeV. 
with the 



the LHC center-of-mass energy we take E, 



76 



We always use the MSTW NLO PDFs 
corresponding value of as{raz) = 0.1202 for the strong 
coupling constant. As our default we use the geometric 
R = \ measure for defining the jets, = 25 GeV for 
the jet veto, and to™' = 200 GeV for the normalization 
range. Our default hard kinematics are = 300 GeV, 
r]j = 0, and Y = 0. Finally, for the scale functions 
Mffi I^b{t)^ 1^j{t), and HSiiT) defined in Sec. IIII Cl the 
central parameter values are given in App. [B) There we 



also discuss the combination of scale variations used for 
estimating the perturbative uncertainties. 



B. Normalization and Convergence 

The unnormalized jet mass spectrum at NNLL with 
our default inputs for the quark and gluon channels are 
shown in Fig. |4(a) As one expects, the gluon jets peak 
at a much higher jet mass than the quark jets. We also 
see that the perturbative uncertainties are quite sizable, 
even at NNLL. 

Normalizing the jet mass spectrum allows one to study 
its shape without contamination from the slow conver- 
gence of the integrated 1-jct cross section, and also re- 
duces the experimental uncertainties significantly. We 
denote the normalized cross section as da/dmj and cal- 
culate it using Eq. p9p where we normalize over the range 
< TOJ < to7'. 



We first study the impact of normalization on the 
perturbative uncertainty. To preserve the normaliza- 
tion, we simultaneously vary the scales in the numera- 
tor and denominator of Eq. (|19p . Comparing the unnor- 
malized cross section at NNLL for the gluon and quark 
channels shown in Fig. 4(a) to the normalized ones in 
Fig. |4(b)] we observe that a substantial portion of the un- 
certainty is related to the integrated cross section rather 
than the shape. In the peak region of the toj spectrum, 
30GcV < Toj < 150 GeV the normalized cross sections 
have a quite reasonable remaining perturbative uncer- 
tainty of ~ 5-10%. 

A big part of the sizable uncertainty in the unnormal- 
ized 1-jet cross section is due to the poor convergence of 
the hard function ion pp ^ H + 1 jet, and thus specific 
to the Higgs process. By keeping the hard function as an 
overall multiplicative factor as in Eq. it cancels ex- 
actly in the normalized cross section for a given partonic 
channel and fixed phase space point (which we use for 
most of our plots) . This cancellation still takes place ap- 
proximately for the integrated cross section summed over 
partonic channels as we show below in Sec. IV Bl Our re- 
sults with fixed kinematics are therefore representative 
of results integrated over the jet phase space. 

The order- by-order convergence of our resummed jet 
mass calculation is displayed in Figs. 4(c) and |4(d)| for 
the gluon and quark jet channels, where results at LL, 
NLL, and NNLL are shown. The various bands overlap 
with those of lower orders, providing direct evidence that 
our scale variations yield a reasonable estimate of the 
higher-order perturbative uncertainties. 

There are several classes of perturbative scale uncer- 
tainties, the "Fixed Order" scale variation that is cor- 
related with the total cross section, the "Beam" scale 
variation from varying /xb and ^Sb ^^at is related to the 
presence of the jet veto, the "Jet" scale variation from 
varying /ij and fisj that is related to the jet mass mea- 
surement, and the uncertainty from "r" that is related 
to the perturbative freedom in the refactorized formula 
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(a) Unnormalized jet mass spectrum for quark and gluon jets 
at NNLL. The uncertainties are sizable even at NNLL. 
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(b) Normalized jet mass spectrum for quark and gluon jets at 

NNLL. Compared to Fig. |4(a)[ the normalization significantly 
reduces the perturbative uncertainties. 
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(c) Convergence of the resummed calculation for gluon jets. 
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(d) Convergence of the resummed calculation for quark jets 
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(e) Individual scale variations that enter the uncertainty estimate 
for gluon jets at NNLL. Shown are the variations relative to 
the central NNLL curve. 
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(f) Individual scale variations that enter the uncertainty estimate 
for quark jets at NNLL. Shown are the variations relative to 
the central NNLL curve. 



FIG. 4: Perturbative uncertainties and convergence for the jet mass spectrum in gg — > Hg and gq — )■ Hq with default inputs. 



for the soft function. For the NNLL resul ts, th ese indi- 
vidual scale variations are shown in Figs. 4(e) and |4(f)| 
for gluon and quark jets respectively. For simplicity we 
combined the uncertainty from varying the jet scale /xj 



and the scale of the jet part of the soft function fj,Sj by 
taking the envelope, and similarly for the beams. It is not 
too surprising that the uncertainties associated with the 
hard and beam scale variations are smaller, since they 
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FIG. 5: Effect of using different jet veto cuts on the jet mass spectrum for gg — >■ Hg. While the unnormalized spectrum on 
the left is directly sensitive to the jet veto cut, this dependence almost completely cancels in the normalized spectrum on the 
right. The same is true for the quark channel, gq — > Hq, and the sum over all partonic channels. 



are mostly common to the numerator and denominator 
of the normaUzed spectrum in Eq. (fT9|) . To obtam the 
total perturbative uncertainty we take the envelope of 
"Jet", "Beam" and "r" uncertainties and combine it in 
quadrature with the "Fixed Order" uncertainty. The to- 
tal uncertainty in the jet mass spectrum is dominated by 
that of the jet and by the soft function refactorization. 



C. Jet Veto and Non-Global Logarithms 

Next we discuss the effect of the jet veto on the jet mass 
spectrum. Our veto is imposed through the variable T'^"', 
rather than a more traditional p^}, since this simplifies 
the treatment of scales in the problem, and allows us to 
make use of a known factorization theorem. This jet veto 
restricts the initial and final-state collinear radiation as 
well as soft radiation. It turns out that the normalized 
jet mass spectrum is fairly insensitive to the value of the 
jet veto cut. 

We start by considering the effect of the jet veto on the 
unnormalized jet mass spectrum, as shown for gg — >■ Hg 
in the left panel of Fig. [5l Decreasing 7"'^^* imposes a 
stronger restriction on the initial-state radiation and re- 
duces the unnormalized cross section. (This reduction 
is less strong for gq — >■ Hq, because quarks radiate less 
than gluons.) As the right panel of Fig. [5] shows, the nor- 
malization removes the majority of the 7"'="* dependence. 
Note that without the refactorization of the soft function 
(see Sec. IIIIBp this cancellation would be spoiled by un- 
physical logarithms. This strong cancellation is also the 
case for the other partonic channels, as well as for their 
sum in — > -I- 1 jet. This insensitivity to T"'^"* also 
remains valid after integrating over the jet phase space, 
as we show below in Fig. [8l We have also studied the 
dependence on T'^"* as well as a standard p™} jet veto 
with Pythia, where we also find a similar insensitivity 
of the normalized jet mass spectrum to the details of the 



used jet- veto variable and cut values. 

Next we turn to our analysis of NGLs, both in the 
unnormalized and normalized jet mass spectra. As ex- 
plained in Sec. IIII Dl we test for the size of the NGLs 
by comparing the cross section with and without these 
terms. The leading NGL is included in fixed-order per- 
turbation theory, on top of which we sum an infinite series 
of global logarithms through the factorization formula. 

In the left panel of Fig. [6] we show the unnormalized 
spectrum for various T*^"' values at NNLL (solid lines) 
and the same spectra including the NGL terms (dotted 
lines). As mentioned earlier, there is a point on the spec- 
trum where the NGLs exactly cancel. This point is at 
mj ~ 50, 110, 165, 300 for T™' = 10, 25, 50, 150 GeV re- 
spectively. For all values of mj shown in this figure the 
effect of the NGL terms is well within the perturbative 
uncertainty [cf. the uncertainty bands shown in Fig. 4(a)| . 

When we normalize the spectrum we arc dividing by 
the cumulant with TOj"*, and the jet- veto dependence 
does not cancel out in the presence of the non-global 
logarithms. There are two types of NGLs in the nor- 
malized result, terms involving ln[mj/(pj.T'^"*)] from 
the numerator and terms involving In 
from the denominator. Therefore for a fixed T"*^"* there 
is no longer a value of mj where all the NGLs will 
vanish. Results for the normalized spectrum with and 
without NGLs are shown in the right panel of Fig. [S] 
The orange band shows the NNLL result without NGLs 
along with its perturbative uncertainty, while the vari- 
ous black lines show the central values for NNLL results 
that have the NGLs included. For the wide range of val- 
ues 25 GcV < r™' < 150 GeV the effect of the NGLs is 
of the same size as the reduced perturbative uncertainty 
in the normalized spectrum. This justifies our assertion 
that the NGLs do not have to be considered as large log- 
arithms for a significant range of cut values, so that our 
NNLL result is complete at this order. In the small mj 
region of the spectrum the resummation of global loga- 
rithms on top of the NGL term provides an appropriate 
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nij [GeV] mj [GeV] 

FIG. 6: Effect of non-global logarithms on the NNLL jet mass spectrum for gg — >■ Hg for different jet veto cuts. Left panel: 
Including the leading NGLs (dashed lines) has a small effect on the unnormalized spectrum, and is well within the perturbative 
uncertainty for a wide range of jet veto cuts. Right panel: The effect of including the leading NGLs (black solid, dashed, and 
dotted curves) on the normalized NNLL spectrum (orange band) is still within the reduced perturbative uncertainty for a wide 
range of jet veto cuts, but start to become important for T"'^"' = 10 GeV. 



Sudakov suppression in the the cross section. For other 
mj values, and 25 GeV < T"'^"' < 150 GeV, the argument 
of the NGL remains between 1/8 and 8, which is the 
range over which we expect that the NGLs do not dom- 
inate over nonlogarithmic corrections, as mentioned in 
the introduction. On the other hand, for T"'^"* = 10 GeV 
one observes that the NGLs become large enough that 
they are no longer contained within the perturbative un- 
certainty, so this value is outside the range of validity of 
our normalized NNLL results (though for this value the 
unnormalized results in the peak region are still valid). 
For this value the argument of the NGL involving m™* 
becomes ~ 13, which is outside of the range mentioned 
above. 

Although we have only explored the gg Hg channel 
at a fixed kinematic point in this section, wc have also 
checked explicitly that the same conclusions about NGLs 
hold when integrating over a kinematic range, and when 
considering quark jets from gq — > Hq. 



V. RESULTS FOR pp ^ H + 1 JET 



In this section we show results for the pp ^ H + 1 jet 
cross section at NNLL, summing the contributions from 
the various partonic channels: gg — > Hg, gq — > Hq, and 
the (small) qq — > Hg. We present results for the depen- 
dence of the jet mass spectrum on the jet kinematics, on 
the choice of jet definition which affects the shape of the 
jets, and on the jet size R. We also compare the mj 
spectrum obtained for a fixed point in the jet kinemat- 
ics to that obtained from integrating over a range of jet 
momenta. 



A. Dependence on Kinematics 

For pp ^ H + 1 jet there are three nontrivial kine- 
matic variables: the transverse momentum of the jet pj., 
rapidity of the jet -qj, and the total rapidity Y of the 
combined Higgs-|-jet system. We show how each of these 
variables affect both the unnormalized and normalized jet 
mass spectrum, which allows us to separate the impact 
of kinematics on the normalization and the shape. 

The falloff of the PDFs at larger x values causes the 
cross section to strongly decrease for increasing p:^ and 
for increasing \r]j\ (for Y = 0). This is shown in Figs. |7(a) 
and |7(c) The dependence on p;^ and ?]j in the corre- 
sponding normalized spectra are shown in Figs. [7(b)] and 
[7(d)| Here wc see that there is a decrease in the height of 
the peak and a compensating increase in the tail height 
as p'^ or \rij\ are increased. Note that for these variables 
there is a marked difference between the total pp ^ H +1 
jet process compared to the individual partonic chan- 
nels (which arc not shown). For each partonic channel 
the peak position of the jet mass spectrum increases as 



m 



peak 



cx 



Pj. and also increases with increasing \r]j\. 
However, at the same time the contribution of gq — >■ Hq 
relative to gg — > Hg is enhanced, and the peak of the 
jet mass spectrum is at lower values for quark jets than 
for gluon jets [see Fig. |4(b)| . These two effects largely 
cancel for pp ^ H + 1 jet, such that the peak position 
is practically unchanged with increasing pj,, whereas for 
increasing rjj a small net increase in the peak position 
remains. 

Note that our ability to calculate the rjj dependence 
implies that it is trivial to impose rapidity cuts in our 
framework. 

The main dependence on the total system rapidity 
Y enters through the shape of the PDFs. causing the 
cross section to strongly decrease with increasing \Y\, as 
Fig. 7(e) shows. (This is also the reason for taking ecu- 
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(d) For pp H + 1 jet the peak position shifts slightly and the 
spectrum slightly broadens with increasing rjj. 
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(f) For pp —)■// + 1 jet the peak position remains stable and the 
spectrum slightly broadens with increasing Y. 



FIG. 7: Dependence on the kinematic variables Pt, ru, and 
for pp ^ H + 1 jet. 



Y for the unnormalized and normalized NNLL jet mass spectra 



tral jets with 7 = for our default value when using 
a single phase space point.) The value of Y also affects 
the shape of the jet mass spectrum, as can be seen in 
Fig. |7(f)[ The jet rapidity relative to the partonic center 



of mass is r] — Y , so one would expect the shape change 
as function of Y to be similar to that as function of 77 j, 
shown in Fig. |7(d)[ The agreement is close but not exact 
because the Y dependence induced by the shape of the 
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FIG. 8: Results for the normalized jet mass spectrum at NNLL for pp H + 1 jet after integrating over 300 GeV < Pt < 
400 GeV, |77j| < 2, and all Y. The left panel compares the spectrum for integrated kinematics (dashed line) to those for fixed 
kinematics with Y = rij — and pj- = 300, 360, 400 GeV (solid lines from top to bottom at the peak). The right panel shows 
the impact of 7"'^"' on the normalized spectrum for integrated kinematics, which is the analog of the comparison in the right 
panel of Fig. [S] for fixed kinematics. 



PDFs differs channel by channel, and thus affects their 
relative weight in the sum over channels. 



B. Integrated Kinematics 

So far we have shown the mj spectra for a fixed point 
in pj., 77 J, and Y. We now consider the impact of inte- 
grating the kinematic variables over a bin with \r]j\ < 2, 
300 GeV <p^ < 400 GeV, and any Y. These kinematic 
ranges are realistic experimentally for jets at the LHC. 

In the left panel of Fig. [5] the jet mass spectrum 
for integrated kinematics is shown by a black dashed 
line, and is compared to three spectra with fixed kine- 
matics shown by solid lines (with V = 77,7 = and 
Pj. = 300, 360, 400 GeV from top to bottom at the peak 
of the spectrum). One observes that the nij spectrum in 
the integrated bin is very close to the mj spectrum with 
y = 77,/ = and near the center of the pj. bin. Thus our 
conclusions made from studies of a single kinematic point 
directly carry over to the results obtained by integrating 
over a phase space bin. 

The one situation where this is not immediately ob- 
vious is the dependence of the normalized cross section 
on the jet-veto cut, T'^"*, shown for fixed kinematics in 
Fig. [5l When we integrate over the kinematic bin the 
hard function, including its Sudakov form factor depend- 
ing on T'^"*, no longer exactly cancels between the nu- 
merator and denominator. Nevertheless, comparing the 
spectra for integrated kinematics and different values of 
T''"', shown in Fig. [HI we sec that the normalized spec- 
trum is still very insensitive to the details of the jet veto 
also after summing over partonic channels and integrat- 
ing over a range of kinematics. (We have also confirmed 
that upon phase space integration the size of the NGL 
effect remains the same as shown in Fig. [S]) 



C. Jet Definitions and Jet Area 

In Sec. |TT] we discussed the various 7V-jettiness mea- 
sures (defined by the Q,) and illustrated the correspond- 
ing size and shape of the jet regions for the geometric 
cases. An illustration of the more irregular regions that 
appear for the invariant mass measure can be found in 
Ref. [131 . We now study how the jet mass spectrum is af- 
fected by these different jet definitions as well as by their 
jet area {R dependence). We start by noting that in the 
iV-jettiness factorization only the soft function is sensi- 
tive to the boundaries of the jet regions. Since the soft 
function only probes the boundaries starting at NLO, the 
first dependence on the jet definition appears at NNLL. 
The nontrivial Qi dependence at this order in the singu- 
lar terms is formally enhanced over the dependence on 
the jet algorithm and jet area in the power-suppressed 
nonsingular terms that are not part of Eq. ([^5]) . 

In Fig. [9] we compare the invariant mass, geometric 
Pt, and geometric E measures (using p = 0.5 for the 
latter two) for three different kinematic configurations 
with Tjj = {0,1,2}, V = 0, and p-^ ^ 300 GcV. For 
?7j = and rjj = 1 the dependence of the jet mass on 
the jet definition is fairly mild (for jets of similar area), 
with the largest visible effect at small mj. As we saw in 
Fig. 1(a) the jet regions of these two geometric measures 
agree exactly for r]j = and differ progressively for more 
forward jet rapidities. This is reflected in the jet mass 
spectrum where the largest differences are observed for 
77,7 — 2. The two bumps in the jet mass spectrum for 
geometric pt at 77,7 = 2 are caused by an increased sep- 
aration between the peaks of the individual quark and 
gluon channels [see Fig. |4(b)"| . 

In Fig. (TU] we show the jet mass spectrum for the ge- 
ometric R measure for various values of the jet radius 
R. A larger jet radius translates into a peak at slightly 
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FIG. 9: Dependence of tlie NNLL jet mass spectrum for pp ^ H + 1 jet on the Ai'-jettiness measure used to define the jets. 
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FIG. 10: Dependence of the NNLL jet mass spectrum for 
the geometric R measure on the jet radius R. On[y the R 
dependence from singular terms in the factorization formula 
is shown here. 



VI. IVIONTE CARLO COIVIPARISONS 

In this section wc study various aspects of the jet mass 
spectrum in Pythia. Although formally the perturba- 
tive accuracy of Pythia is significantly lower than that 
of our NNLL calculation, it is also well known that after 
sufficient tuning Pythia is able to reproduce the shape of 
many jet observables. Here we are particularly interested 
in testing the impact on the jet mass spectrum from using 
different hard processes, using different jet algorithms, 
and from adding hadronization and underlying event (the 
latter being described by Pythia's multi-parton interac- 
tion model). We also perform a comparison between our 
calculation and Pythia for the same geometric R = 1 
A^-jcttiness jets used in our analysis. Finally we compare 
our exclusive 1-jet mj calculation with the inclusive jet 
mass spectrum measured in pp — > jets by ATLAS. Wc al- 
ways use PythiaS with its default tune 5 ("Tune 4C"), 
which as we will see provides a good description of the 
ATLAS jet mass data. 



A. Hard Process and Jet Algorithm Dependence in 
PYTHIA 



higher mass and slightly larger tail. Note that the biggest 
difference between the jet mass spectra for different val- 
ues of R will be at large mj in the tail of the spectrum, 
since the size of the jet puts an effective upper bound- 
ary on its mass nij ^ pj^R/V2. ^ At this boundary the 
jet mass spectrum has to fall off rapidly. This bound- 
ary is seen in Pythia and LHC data but is not included 
in our calculation, because we have not yet incorporated 
the nonsingular contributions to the cross section which 
are important for accounting for this effect and making 
realistic predictions in this part of the tail of the distri- 
bution. 



^ For a jet of fixed three momentum centered at 77 = 0, tlie abso- 
lute bound on the jet mass is < cos^ R — I), which 
is reached vyhen the jet consists of two energetic narrow sub- 
jets located at (??, fli) = (0, ±R). Assuming a uniform energy 
distribution in {rj, 0)-space, leads to a more practical bound 
m2 < p-pR2/2. 



We start by investigating to what extent the jet mass 
spectrum depends on the underlying hard process in 
Pythia. In Fig. [TT] we show the spectrum for a gluon 
jet from gg — > gg (solid) and from gg — )- Hg (dotted), 
demonstrating that in Pythia there is essentially negligi- 
ble process dependence for individual partonic channels. 
This is true both at the partonic level (blue curves with 
peak on the left) and after including hadronization and 
multiple interactions (red curves with peak on the right). 
In reality one expects some differences from the hard pro- 
cess due to the additional soft radiation produced with 
more available colored particles, and from the different 
color flow, where in particular gg — )■ gg involves a matrix 
of color channels with nontrivial interference. These ef- 
fects may not be sufficiently described by Pythia so one 
should not conclude that the hard process dependence on 
the jet mass spectrum is as small as is shown. 

Next, we look at the difference in Pythia between 
the jet mass for exclusive and inclusive jet production. 
We use the process gg — Hg, imposing the jet veto 
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FIG. 11: The gluon jet mass spectrum in Pythia does not de- 
pend on the underlying hard process producing the jets. This 
is true both for partons (left peaks) and with hadronization 
and underlying event (right peaks). 
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FIG. 13: Comparison of the anti-fcT, CA, and geometric R 
jet algorithms in Pythia. 



different jet measures for 1-jettiness jets. The difference 
between R = \ and R = 1.2 for anti-k^ is larger than 
that observed in our calculation using geometric R jets 
in Fig. \W\ In Pythia the difference between R ~ I 
and R = 1.2 becomes smaller when 7"="' is decreased, 
since with a stronger jet veto less additional radiation is 
present that would be absorbed by larger jets. To be spe- 
cific, the 15% difference in the peak heights for anti-k^ 
with i? = 1 and i? = 1.2 for T™* = 25 GeV reduces to 7% 
for 7""="* = 5 GeV. From these results we conclude that 
in our NNLL calculation, the R dependence of the non- 
singular terms that are not included in our analysis may 
well be numerically comparable or larger in size than the 
formally leading R dependence that we have computed. 



FIG. 12: Comparison of the normalized jet mass spectra for 
exclusive and inclusive jet samples in Pythia. 



B. Comparison of NNLL with PYTHIA 



j-cut _ 25 GeV to obtain two exclusive samples, and 
using no jet veto for our inclusive sample. The resulting 
normalized jet mass spectra are shown in Fig. 1121 The 
difference between 7"^="* = 25 GeV (our default value) 
and the inclusive case is small, allowing our calculation 
to be compared to inclusive spectra. The difference is 
slightly larger for 7"'^"* = 10 GeV and increases signifi- 
cantly for smaller values of T^"*. However, we will not 
consider such strong jet vetos, as they lead to large NGLs 
(see SecHVC]). 

In Fig. [13] we compare the jet mass spectrum from 
Pythia for different jet algorithms, specifically our 1- 
jcttiness R = 1-algorithm, Cambridge-Aachen with R ~ 
1, and anti-kT with i? = 1 and i? = 1.2 (t^I- To stay 
close to a calculation for a single phase space point, we 
restrict the jet to a narrow pr and rapidity bin, and im- 
pose a veto using 7""^"' = 25 GeV. The differences be- 
tween the R = 1 curves are within the size of the uncer- 
tainty band from our NNLL calculation in the same phase 
space bin. This result agrees with the small differences 
observed in each of the panels of Fig. [Hlfrom comparing 



A comparison between our NNLL calculation and par- 
tonic Pythia results for gg — > Hg are shown in Fig. [T4l 
both using the geometric R = 1 jet definition. The peak 
positions in both cases agree very well. To ensure that 
this is not an accident and that the peak position in 
Pythia does not depend on the PDF set used by our de- 
fault tune, we checked that an alternative tune (number 
10, which is based on our default Pythia tune but uses 
MSTW2008 LO PDFs) only shifts the peak by a smaU 
amount, similar to the small difference in peak positions 
between Pythia and our NNLL calculation. However, as 
seen in Fig. [141 the NNLL calculation has a lower peak 
and a correspondingly higher tail. Since the spectrum is 
normalized these two effects are related, namely higher 
values in the tail must be compensated by a lower peak. 
There are several possibilities that may account for this 
difference. Due to the stability of our order-by-order re- 
sults in Fig. 4(c) it is unlikely to be related to the lower 
order accuracy of Pythia's LL parton shower resumma- 
tion. Most likely the differences are related to the fact 
that we have not yet included nonsingular contributions 
to the spectrum which are important in the tail region. 
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FIG. 14: Comparison between our NNLL calculation and par- 
tonic Pythia for the gg — ^ Hg channel. Both results use 
geometric R = 1 jets and the same kinematic cuts. 



in particular for the spectrum to fall off rapidly enough. 
Due to the fact that the results are normalized, this mis- 
match in the tail then also leads to a disagreement of the 
peak heights. Thus we expect that the inclusion of the 
nonsingular contributions will reduce this difference. 



FIG. 15: The nonperturbative hadronization correction in 
Pythia is well described by a shift in m^. 



in Pythia, by comparing the results with and without 
hadronization. As demonstrated in Fig. 1151 a shift with 
the choice 17 = 0.8 GeV works very well, in reasonable 
agreement with the fi = 1.0 GeV found earlier in Ref. [l3| 
for the inclusive > 1 jet cross section. 



C. Hadronization in PYTHIA 

We now explore the effect of hadronization on the jet 
mass spectrum using Pythia. In the factorization for- 
mula the hadronization is encoded through nonperturba- 
tive corrections in the soft function 5 at a scale ^ Aqcd , 
which must be separated from perturbativc corrections 
at the soft scale fis ^ itt-'j/p't- For e+e^ 2 jets there 
is an analytic understanding of the analogous nonpertur- 
bative corrections originating in Refs. j78l - l8ll | as well as 
a modern understanding in terms of field theory opera- 
tors [TqI. l82| - [8^ . For these processes, as soon as the rel- 
evant soft scale fxs is perturbative, the nonperturbative 
corrections can be power expanded in Aqcd//^s, and the 
dominant power correction simply shifts the event shape 
distribution, e — > e — fle/Q. In the case at hand, the 
nonperturbative soft function is built from more than 
two Wilson lines, so the description of the power cor- 
rections becomes iTiore complicated. Nevertheless, for 
a given kinematic configuration we still expect that the 
dominant effect will be described by a shift involving a 
parameter il ^ Aqcd- For a jet mass mj ~ PjPj this 
shift occurs due to nonperturbative soft radiation caus- 
ing a shift in the small momentum pj", so it takes the 
form 

m^j ^ m^j -2np^j,R. (35) 

The factor of R accounts for the fact that there is a de- 
creased amount of soft momentum contamination in the 
jet for decreasing R [sH]. It is straightforward to test 
whether this shift agrees with the hadronization model 



D. Underlying Event and ATLAS Data 

In Pythia the effect of the underlying event is mod- 
eled by multiple partonic interactions, and its effect on 
the jet mass spectrum is more proirounccd than that 
of hadronization. This is shown in Fig. [16] where we 
plot the jet mass spectrum for inclusive pp — >■ jets from 
Pythia at parton level, including hadronization, and in- 
cluding hadronization and multiple interactions. Also 
shown are the corresponding ATLAS data from Ref. [2^ , 
where the uncertainty bars are from linearly combining 
the statistical and systematic uncertainties. This chan- 
nel is dominated by the copious pp — > dijet production 
at the LHC. We use the same inputs and cuts as AT- 
LAS, namely i?cm = 7TeV, anti-k^ jets with i? = 1, 
|?7j| < 2, and consider both 300 GeV < p^ < 400 GeV 
and 500 GeV < < 600 GeV. The shift to the peak lo- 
cation from hadronization is of similar magnitude as that 
for gg — Hg in Fig. [151 namely ~ 3.0 GeV for gg — ?> Hg 
compared to ~ 8.0 GeV for the 300 GeV < < 400 GeV 
inclusive jets which have a slightly larger average p^. 
For the inclusive pp — >■ jets in Pythia the additional 
shift to the peak location from the underlying event is 
~ 17.4 GeV. The final Pythia results agree weU with 
the ATLAS data for both pj. bins. In a NNLL calcula- 
tion the effect of hadronization and part of the effect of 
the uirderlying event will be captured by corrections to 
the soft fuirction, but it is not clear if hadronic correc- 
tions in the multi-jet soft function will fully capture the 
effect of the underlying event. 

Given that Pythia agrees well with the ATLAS inclu- 
sive dijet spectrum, one might wonder what the purpose 
of a higher-order NNLL dijet calculation would be. An 
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FIG. 16: Comparison of the Pythia jet mass spectrum for inclusive pp — >■ jets to the corresponding ATLAS data [26|. 
Pythia results are shown at parton level (dotted), including hadronization (dashed), and including hadronization and multiple 
interactions (solid). The final Pythia results reproduce the data well. 
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FIG. 17: Comparison of our exclusive NNLL calculation with ATLAS inclusive jet mass data [2a ]. The peak position of our 
gluon jets from gg — >■ Hg agrees remarkably well with the inclusive dijet data. For the ATLAS date there is presumably a 
shift to lower values due to quark jets which is compensated by a shift to higher values due to hadronization and multiple 
interactions. 



advantage of our calculational framework over Pythia is 
that it follows from first principles and docs not involve 
the modeling and tuning present in Pythia. Specifically, 
the input to our calculation is limited to a^imz), the par- 
ton distributions functions, and simple soft function pa- 
rameters like n for the hadronic effects. Furthermore, we 
have a rigorous estimate of the higher-order pcrturbative 
uncertainty from scale variation, as well as from order- 
by-order convergence, which enable us to fully asses the 
reliability of the result. Finally, it should be emphasized 
that our calculation is fully analytic (up to the numer- 
ical convolution with the PDFs) and hence provides an 
analytic QCD calculation of an LHC spectrum for jets. 

To the extent that the normalized jet mass spectrum is 
independent of the hard process and independent of using 
an inclusive or exclusive jet sample, which Pythia seems 
to suggest in Figs. [TT] and [T^ a comparison between jet 
mass spectra involving different hard processes and with 
and without jet veto cuts is appropriate. The approx- 



imate hard process independence only holds separately 
for gluon or quark jets, which themselves have fairly dif- 
ferent jet mass spectra, see Fig. |4(b)[ Therefore when 
varying the hard process we expect the dominant change 
in the jet mass spectrum to be related to the process 
dependent fraction of quark and gluon jets produced. 

In Fig.[T7]we compare our NNLL result for — > H + 1 
jet and for gg — >• Hg to the ATLAS data for pp jets. 
Recall that the peak location of the NNLL H+1 jet calcu- 
lation matches well with that from Pythia, see Fig. [141 
Because of the significant contribution from quark jets 
the -ff + 1 jet spectrum peaks to the left of the spectrum 
from dijcts. On the other hand, the peak location with 
pure gluon jets {gg Hg) agrees remarkably well with 
the data on dijets. From the results already obtained 
above, we expect only small differences (comparable to 
the ATLAS error bars) for effects related to the choice 
of the jet algorithm, the choice of inclusive versus exclu- 
sive jets, or the choice of looking at gluon jets in dijcts 
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or in Higgs production. On the other hand there will be 
a more significant shift of the spectrum to the left from 
quark channels in the dijet production, and a shift to the 
right from adding hadronization and underlying event, 
neither of which is included in the solid red curve. The 
agreement between peak locations seems to indicate that 
these two effects largely compensate for one another. Fi- 
nally, there will be an effect related to the fact that there 
are nontrivial color correlations in gg gg which are 
not present in gg — >■ Hg (these effects are not apparent 
in Pythia, see Fig. [TT|l . 

One may also look at the peak heights in Fig. [171 for 
which the agreement is not as good. As described earlier, 
this effect is related to the fact that we have not yet 
included nonsingular corrections. These corrections are 
known to decrease the tail to enable it to rapidly fall off 
by rrij ^ and they also affect the peak directly 

through the normalization. Since with additional work 
these can be included in future results the difference in 
peak heights is not of too much concern. 

Finally one may also compare the results in Figs. 
and [17] for the 300 GeV < < 400 GeV and 500 GeV < 
Pj. < 600 GeV bins. For dijets the peak location moves 
to higher mj with increased p;^, unlike for pp H + 1 
jet, again indicating that gluon jets likely dominate. The 
conclusions from the comparison with Pythia and the 
contrast to our NNLL calculation remains the same for 
these two ranges of py. 



VII. CONCLUSIONS 

In this paper we calculated the jet mass spectrum for 
pp ^ H + 1 jet to NNLL order. For this exclusive 1-jct 
cross section we veto additional jets with the 1-jcttincss 
event shape, and used the 1-jettiness factorization for- 
mula in terms of hard, beam, jet, and soft functions to 
obtain our results. For the normalized jet mass spec- 
trum the remaining higher-order perturbative uncertain- 
ties from scale variation are at the ~ 5 — 10% level at 
NNLL order, and in addition our results exhibit excel- 
lent order-by-order convergence. 

The normalized NNLL spectrum is quite insensitive to 
the jet veto over a wide range of values, even when ac- 
counting for non-global logarithms. Thus in our frame- 
work non-global logarithms can be accurately treated as 
fixed-order contributions to the soft function, upon which 
additional global logarithms are automatically resummcd 
in the factorization framework. An essential ingredient in 
the resummation of the global logarithms was the rcfac- 
torization of the soft function, which we demonstrate is 
required to avoid introducing spurious leading logarithms 
in certain regions of phase space. Our treatment of the 
NNLL exclusive cross section with a jet veto has signif- 
icantly smaller non-global logarithmic terms when com- 
pared to the size of these terms observed in the earlier 
inclusive NLL analysis in Ref . [l3l , and the earlier inclu- 
sive partial NNLL analysis in Ref. [3l|. Finally, we note 



that in Pythia the inclusive jet mass spectrum and the 
exclusive jet mass spectrum with our default jet veto are 
essentially identical. 

Utilizing our calculation we investigated the depen- 
dence of the jet mass spectrum on various parameters 
of the exclusive jet cross section. Part of the power of 
our framework is that the factorization formula is fully 
differential in the jet kinematics (pj., rjj, and Y), allow- 
ing us to vary the definition of the jets and the jet area, 
and can be easily separated into quark jet and gluon jet 
channels. As expected we find that the spectrum peaks 
at larger mj values for gluon jets than for quark jets. For 
a given partonic channel the factorization framework pre- 
dicts little sensitivity to the underlying hard process, and 
this result is also found to be the case in Pythia. The 
main process dependence is therefore the relative mix of 
quark and gluon jets. The peak of our NNLL mj spec- 
trum moves to the right for larger p:^ and for larger jryjl, 
but more so for the individual partonic channels than for 
pp — > H+1 jet, where the change to the mix of quarks 
and gluons provides a compensating effect. The complete 
description of the various kinematic variables also makes 
it trivial to implement rapidity cuts. For a bin \rij\ < 2 
and a not too large bin in p j. , we find that the integrated 
NNLL result is very consistent with the NNLL result for 
fixed kinematic variables taken at the center of the bin. 

Varying the jet definition with fixed jet area leads to 
very little change in the jet mass spectrum, both for var- 
ious jet definitions in our NNLL result and for anti- /ct, 
CA, and geometric- i? jets in Pythia. This suggests that 
there are only small differences between the spectrum 
for 1-jettiness jets and traditional jet algorithms. On the 
other hand, Pythia exhibits a larger dependence on the 
jet radius R than our NNLL results. This is presumably 
due to the fact that we have not included nonsingular 
terms in our analysis. These terms become important in 
the tail region and are also needed to ensure that the jet 
mass spectrum dies off at mf^^ ^ pj,R/y/2. This leads 
to a larger tail in our NNLL spectrum than in Pythia, 
and correspondingly a smaller peak height in the nor- 
malized NNLL result. On the other hand, the peak lo- 
cation agrees very well between our NNLL calculation 
and Pythia. An analysis of these additional nonsingu- 
lar terms will be carried out in the future. 

We investigated the dependence of the jet mass 
spectrum on hadronization and underlying event using 
Pythia. Hadronization is very well described by a 
shift to the mass spectrum, TOj m ^ — {2Rp^)Vl with 

^ Aqcd, which is the anticipated result from non- 
perturbative soft gluon contributions in our factorization 
formula's soft function. In Pythia the underlying event 
is modeled by multiple partonic interactions and leads 
to a somewhat larger shift to the spectrum than for 
hadronization. It plays an important role in obtaining 
agreement with the ATLAS jet mass results for inclusive 
dijets. Comparing our results to ATLAS we find that 
the NNLL pp — >■ H+\ jet spectrum peaks to the left of 
the dijet data, whereas the NNLL gg — >■ Hg spectrum 
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peaks in the same location. The comparison made so far 
with the ATLAS data is promising. The extension of our 
NNLL calculation to pp — !• dijets is completely feasible 
using 2-jettiness, and it will be interesting to see to what 
extent the contributions from quark channels, color mix- 
ing, and hadronization and underlying event will affect 
this comparison with the data. Theoretically, the only re- 
maining challenge to a complete comparison appears to 
be incorporating the effect of the underlying event from 
first principles rather than relying on its modeling via 
Monte Carlo. 



Appendix A: Perturbative Inputs 

In this section we collect the fixed-order ingredients 
and evolution kernels for evaluating the jet mass cross 
section for pp ^ H + Ij in Eqs. ^ and (HH) at NNLL 
order. We first give expressions for the hard, jet, beam 
and soft functions at next-to-leading order. This is fol- 
lowed by the evolution kernels and the coefficients that 
they depend on. 
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1. Hard Function 

The hard functions ff^ for the various partonic chan- 
nels K that contribute to pp ^ H +1 jet can be obtained 
from the finite part of the helicity amplitudes A deter- 
mined in Ref. 1561 . 



Hggg{{q'l},^^H) = 
Hqgqi{qt},f^H) = 



IQasiiiHYCiCF 



97rw2 [2(iV2 - 1)]2 
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Hqqgi{q'^},f^H) = ^ggg ({gf} , ) • 



(Al) 



The factors of l/(2iVc) and 1/[2(A^2 _ i)] arise from averaging over the spins and colors of the colliding quarks and 
gluons. The arguments of a helicity amplitude A have the form iJ*, where i denotes the momentum , t denotes 
the parton type, and h denotes the helicity of this particle. Only in the helicity amplitudes will we use an outgoing 
convention for all these quantities, to make crossing symmetry direct. This implies that if we want to convert to the 
convention used in the main text, then the Sij's in the helicity amplitudes below will pick up additional minus signs 
if one of the particles i and j is in and the other is out. The amplitudes that enter in Eq. (jAip are given by 
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Here we use the shorthand notation 
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The Cn{x) integrate to zero if the range in x is [0, 1]. 



2. Jet Functions 



The one-loop jet functions are given by |58l-l6i 



Jq{s,^ij) = 5{s) + 



ois{iJ'j)Cf r 2 / s 



27r 



2^l( 2 
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+ [(^- y)Ca + ^/3o]'5(s)}, (A4) 
where the plus distributions £„ are defined as 
r0(x)ln"x 



3. Beam Functions 

The beam functions can be expressed in terms of stan- 
dard gluon and ouark PDFs using an operator product 
expansion [61I, |8g| , 
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The one-loop matching coefficients are [62, ISJ] are 
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The splitting functions in this equation are defined as 

Pqq{z) = Coil - Z){1 + Z^) , 



Pgg{z)^eil-Z)[{l-Zr + Z'], 
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Pgg{z)^2CQ{l- z)z + 2e{l- z) 



+ z{l~z) 



Pg,{z)^e{\-z) 



l + {l-zf 



(A8) 



freedom in this refactorization, and that the correspond- 
ing uncertainty is probed by varying the parameter r in 
Eq. dSIl). 

Up to NLO the 1-jettiness soft function is given by 



4. Factorized Soft Function 

We now give expressions for the TV-jettiness soft func- 
tion, showing expUcitly how the factorization in Eq. (|24p 
is implemented. We remind the reader that their is some 
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From the NLO calculation in Rcf. [47| we obtain 
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Following the discussion in Sec. IIII B( we give expres- 
sions for the factorized evolution of the hard function, 

ii,{X(i;},{y.,}) = R^{{q>t},m) n UH^iX^}.m.P'^). 



n( 
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VHil^HyfJ'i) ^ = . (A13) 
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Here the products over i and j run over all colored par- 
ticles, with corresponding flavor Hi and Hj. For each 
channel contributing to — > H + Ij there is only 
one color structure so • Tj is simply a number [see 
Eq. (|A12p ]. The functions Kr, i]r and are given be- 
low in Eq. ([XTH)) . 

The solution of the RG evolution of the jet function is 
given by m m 113, HI 



U.j^^{s,fij,f^) , 

r(l-hr?» 

The plus distribution is defined as 
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General relations for the rescaling and convolutions of 
Cn{x) in Eq. (jASp and can be found in App. B 
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of Ref. [7^ . The renormalization group evolution of the 
beam functions is identical [62j and can be obtained from 
the above expressions by replacing Ji{s, fi) — ?> Bi{t, x, fi). 
We do not give the evolution of the soft function, as it 
is not needed for evaluating Eq. ([28|) . It can be obtained 



from the evolution of the hard function and beam func- 
tion by using the /^-independence of the cross section. 

The functions A'r(/io,M), f?r(/^0:M)i K^(fio,^) in the 
above RGE solutions at NNLL are given by, 
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Here, r = as{fi)/as{fJ.o) and the running coupling at the 
scale fi is given in terms of that at the reference scale fio 
by the three-loop expression 
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where X = l-f as(Mo)/3o ln(At/A'o)/(27r). 

6. RGE Coefficients 

Up to three loops, the coefficients of the beta func- 
tion [H, [9^ and cusp anomalous dimension [6^ |6^ in 
MS are 
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Up to two loops, the MS non-cusp anomalous dimen- 
sion for the hard function [qiI,!!^ and jet and beam func- 
tions ^ E2, 113, [t^ are 
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Appendix B: Running Scales 

We now present the remaining ingredients that enter 
in the running scales in Sec. IIII CI First of all, /ij-un is 
defined as 
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The expressions for a and h follow from demanding that 
A*run(T) is continuous and has a continuous derivative. 
For our central scale choice we use 

^ = Q , = es, = , /io = 2GeV , 
5GeV 



Q 



T2 = 0A, Ts^O.e, r = 0.2. (B2) 



To estimate the perturbative uncertainty we vary the 
above parameters within reasonable ranges. Since the 
cross section is most sensitive to fi, e^, eg^ and r, we 
restrict ourselves to the following separate variations, 
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Following our discussion in Refs. [93|,|9J], we take the en- 
velope of variations b) through f ) and add this in quadra- 
ture with variation a). 
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